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ABSTRACT 

We follow up our deep ROSAT imaging survey of the Pleiades (Stauffer et at. 1994) with an analysis of the 
spectral and temporal characteristics of the X-ray-bright stars in the Pleiades. Raymond & Smith (1977) one- 
and two-temperature models have been used to fit the position-sensitive proportional counter (PSPC) pulse- 
height spectra of the dozen or so brightest sources associated with late-type Pleiades members. The best-fit 
temperatures suggest hot coronal temperatures for K, M, and rapidly rotating G stars, and cooler tem- 
peratures for F and slowly rotating G stars. In order to probe the many less X-ray-luminous stars, we have 
generated composite spectra by combining net counts from all Pleiades members according to spectral type 
and rotational velocity. Model fits to the composite spectra confirm the trend seen in the individual spectral 
fits. Particularly interesting is the apparent dependence of coronal temperature on L^jL^ A hardness-ratio 
analysis also confirms some of these trends. The PSPC data have also revealed a dozen or so strong X-ray 
flares with peak X-ray luminosities in excess of ~ 10 30 ergs s _i . We have modeled the brightest of these flares 
with a simple quasi-static cooling loop model. The peak temperature and emission measure and the inferred 
electron density and plasma volume suggest a very large scale flaring event. The PSPC data were collected 
over a period of - 18 months, allowing us to search for source variability on timescales ranging from less than 
a day (in the case of flares) to more than a year between individual exposures. On approximately year-long 
timescales, roughly 25% of the late-type stars are variable. Since the Pleiades was also intensively monitored 
by the imaging instruments on the Einstein Observatory , we have examined X-ray luminosity variations on the 
10 yr timescale between Einstein and ROSAT and find that up to 40% of the late-type stars are X-ray vari- 
able. Since there is only marginal evidence for increased variability on decade-long timescales, the variability 
observed on long and short timescales may have a common physical origin. 

Subject headings: open clusters and associations: individual (Pleiades) — stars: coronae — stars: late-type — 
stars: rotation — X-rays: stars 


1. INTRODUCTION 

The Pleiades is an important laboratory for the study of 
main-sequence stellar evolution because of its youth and prox- 
imity. With an estimated cluster age ~ 70 million years, many 
of the intermediate- and low-mass dwarfs are still rapidly 
rotating, providing clues to the evolution of angular momen- 
tum and lithium depletion on the main sequence. Rapid rota- 
tors are also expected to be copious X-ray emitters, and, at a 
distance of ~ 127 pc, a large number of intermediate- and low- 
mass cluster members are detectable with X-ray telescopes. 
Consequently, the Pleiades was one of the primary open clus- 
ters observed by the Einstein Observatory (Caillault & Helfand 
1985; Micela et al. 1990) and ROSAT (Stauffer et al. 1994; 
Micela et al. 1995). The analyses presented here represent a 
continuation of the Stauffer et al. (1994, hereafter Paper I) 
work. In that study, more than 300 X-ray sources were 
detected, of which 171 are associated with certain or probable 
cluster members. We follow up that spatial analysis with spec- 
tral and temporal analyses of X-ray sources associated with the 
late-type Pleiades stars. 

The study of stellar X-ray spectra has been difficult because 
of the relatively poor energy resolution of the Einstein imaging 
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proportional counter (IPC) and the ROSAT position-sensitive 
proportional counter (PSPC). Nevertheless, many insights 
have been gained. For example, Schmitt et al. (1990), in a com- 
prehensive study of IPC pulse-height spectra of late-type field 
stars, found that one- and two-temperature Raymond-Smith 
(1977, hereafter RS) optically thin plasma models fitted to the 
observed spectra were largely influenced by the signal-to-noise 
ratio of the data and that continuous emission measure dis- 
tribution models were equally adequate and, perhaps, more 
physically meaningful. Schmitt et al. (1990) also found that M 
dwarfs display evidence of possessing high-temperature gas 
(T > 10 7 K) in addition to lower temperature (T * 3 x 10 6 K) 
material; the F and G dwarfs seemed to have only the low-T 
component. 

In a study of the Hyades stars, Stern et al. (1994) found that 
isothermal coronal plasma models do not fit the observed 
PSPC pulse-height spectra of the 10 brightest sources in the 
cluster (the only ones for which fits were attempted). Fits with 
models having two temperature components suggest that the 
higher of the two temperature components of the M dwarf 
coronae are generally hotter than those of the F-G dwarf 
coronae, similar to the results found by Schmitt et al. (1990). 
Stern et al.’s interpretation of their results is that increased 
volumetric coronal heating, either in a steady fashion or in 
quasi-continuous flaring, may be producing the trend seen 
along the main sequence. They warn, though, that both their 
sample and that of Schmitt et al. are biased toward the high 
end of the X-ray luminosity function. 

Although the stars in the Pleiades are intrinsically more 
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X-ray luminous than those in the Hyades, the fact that they are 
3 times more distant means that their X-ray fluxes are similar 
and, since the ROSAT PSPC exposure times for both clusters 
were approximately equal, the total counts observed for a 
typical individual star are approximately the same. Hence, as 
in the case of the Hyades, only a handful of Pleiades stars have 
sufficient counts with which to perform model fits to their 
X-ray spectra. As in the Hyades, this means we are only 
probing the most X-ray-luminous sources. In order to include 
fainter Pleiades members, we have constructed composite 
X-ray spectra for different spectral types as a function of pro- 
jected rotational velocity (see § 2 below). In this manner, we 
hope to be able to determine in a statistically significant 
manner whether the coronal temperatures are dependent on 
spectral type (as implied by the previous studies) and/or on 
rotation (in some as yet to be determined manner). In addition, 
we can confirm (or refute) the expectation (Schmitt 1994) that 
the temperatures of the model components increase during a 
flare, since we observed one large flare from a star (H II 2147) 
for which its quiescent state is also sufficiently active to acquire 
a statistically valid model fit to its X-ray spectrum. 

Preliminary temporal analyses of Pleiades X-ray data have 
been made by Schmitt et al. (1993), Stauffer et al. (1994), and 
Micela et al. (1995). Schmitt et al. (1993) compared ROSAT 
all-sky survey data with IPC data and suggested that the 
observed source variability in the decade between Einstein and 
ROSAT may be evidence for solar-like coronal activity cycles. 
In Paper I we treated the variability in only a cursory manner, 
determining that, as an ensemble, there is not great variation 
(less than a factor of ~ 4) in the X-ray emission from these stars 
over the decade between observations. Micela et al. (1995) find 
that ~90% of the Pleiades stars are not variable (within a 
factor of ^ 2) on a timescale of approximately 6 months; they 
also find three stars which flared during their PSPC obser- 
vation. 

Hence, in addition to the spectral analysis, we have followed 
statistical techniques described by Schmitt et al. (1993), Gagne 
& Caillault (1994), and Gagne, Caillault, & Stauffer (1995) for 
measuring the variability of X-ray emission over many differ- 
ent timescales, ranging from minutes to more than a decade. 
We compare the medium-term variability (timescales of ~ 1 yr) 
of the Pleiades stars with that of younger stars (e.g., those in 
Orion or p Oph) and older stars (e.g., those in the Hyades), in 
terms of both frequency and strength of variation. By compar- 
ing the ROSAT PSPC and Einstein IPC count rates, we deter- 
mine whether signatures of long-term solar-like dynamo 
activity cycles can be observed among the Pleiades stars. 

2. PSPC SPECTRAL ANALYSIS 

The ROSAT PSPC data for the core of the Pleiades cluster 
used in this paper have been described in detail in Paper I. The 
data consist of three pointed observations with individual 
exposure times ranging from 20.5 to 27.3 ks and covering a 
total of ~4.5 deg 2 . A total of 317 X-ray sources were detected 
in this area, 171 of which have been associated with known or 
probable Pleiades cluster members. 

The ROSAT X-ray telescope and PSPC have been described 
in detail by Pfeffermann et al. (1986) and Triimper (1992), 
respectively. The limited energy resolution (E/AE ~ 2.5 at 1 
keV) of the PSPC has allowed us to construct background- 
subtracted pulse-height spectra of the brighter (net counts 
> 500) sources. In addition, for faint and bright sources associ- 
ated with late-type Pleiades members, we have combined the 


net counts according to spectral type and rotational velocity, 
thereby allowing us to look for differences in the PSPC spectra 
of rapidly and slowly rotating F, G, K, and M stars. 

2. 1 . PSPC Spectra of Individual Cluster Stars 

Of the 171 X-ray sources associated with a single Pleiades 
cluster member, 1 5 were detected with 500 or more net counts 
in at least one of the three PSPC exposures. In the case of the 
bright X-ray flare star H II 2147 (see § 3), sufficient counts have 
been collected in and out of flare to construct both flare and 
quiescent spectra. 

Counts have been extracted from 90% power-radius source 
circles and from concentric background annuli. The local 
background subtraction procedure we used is described in 
detail in Paper I. In particular, the background regions have 
been masked from nearby sources and the effective exposure 
times for each source have been estimated from broadband 
PSPC exposure maps. 

In some cases a source has been detected in two exposures. 
In order to improve the signal-to-noise ratio of the X-ray 
spectra of these stars, the counts in both the source and the 
background regions from the two exposures have been co- 
added and net rms errors have been determined. Since the 
energy response of the PSPC depends on the source’s position 
on the detector, the net energy response of the detector for the 
co-added spectrum was determined by taking a weighted mean 
of the individual responses. 

The resulting PSPC pulse-height spectra were fitted with 
one- and two-temperature RS models within PROS. (The 
PROS software is provided by the Harvard-Smithsonian 
Center for Astrophysics.) The hydrogen column density was 
fixed at 2 x 10 20 cm~ 2 , consistent with the mean color excess 
E(B— V) ~~ 0.04 found in the core of the cluster (see Paper I 
and references therein). When allowed to float, the best-fit 
column densities were always in the range (1. 6-2.5) x 10 20 
cm' 2 , thereby supporting our assumption of approximately 
constant N H , at least for the purpose of fitting PSPC pulse- 
height spectra. In all cases, we find acceptable fits assuming 
cosmic abundances. That is not to say, however, that these 
coronae have cosmic abundances. The energy resolution and 
signal-to-noise ratios of the PSPC pulse-height spectra of the 
Pleiades stars are simply too poor to constrain coronal abun- 
dance values. 

As was found by Schmitt et al. (1990) for stars detected with 
the Einstein IPC, we find that single-temperature models ade- 
quately fit the observed spectra when the total counts in the 
spectrum are low. For example, for five of six stars with fewer 
than 1000 net counts, a single-temperature model provides an 
acceptable fit, i.e., xl ^ 1* This, of course, is an artifact of the x 2 
statistic and indicates that the low-count spectra cannot con- 
strain multitemperature models because of the large relative 
errors in each spectral bin. For all sources with net counts 
above 1000, two-temperature RS models were required. 

In Table 1 we present the results of the RS model fits to the 
15 sources with >500 net counts associated with a known 
Pleiades cluster member. In column (1) we list the X-ray source 
number from Paper I; in columns (2)— (4) we list the name, 
(B — V) 0 color, and rotational velocity of the Pleiades stars. 
The remaining columns list the results of the spectral fits. In 
columns (5)— (10) we list the total number of net counts in the 
0. 1-2.0 keV band, the best-fit values of x 2 , the volume emission 
measure of the first temperature component (EM,), the tem- 
perature of the first component (fcT,), EM 2 , and kT 2 , respec- 
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tively. Stars with no data in columns (9) and (10) were well 
fitted with a single-temperature model. We note that the FSPC 
pulse-height spectra have been corrected for vignetting, the 
off-axis decrease in telescope sensitivity. 

The parameters in Table 1 are not meant to indicate abso- 
lute coronal temperature distributions. We use a two- 
component temperature distribution simply because these 
PSPC data cannot constrain more sophisticated models. In 
particular, continuous emission measure distributions do not 
provide better fits to the data. The best-fit temperatures and 
volume emission measures in Table 1 should, however, allow 
us to compare the relative distributions of coronal plasma 
among the different stars observed by ROSAT . As expected, 
during the H II 2147 flare we observe higher temperatures and 
much larger volume emission measures than during quies- 
cence. 

In Figure 1 we plot kT 2 against kT { for the 10 stars in Table 
1 that have been fitted with the two-temperature model, 
including H II 2147 in and out of flare. The X-ray flare stars 
noted in Table 1 are indicated with a cross. The rapid rotators 
(y sin i > 16 km s' 1 ) are indicated with filled symbols, while 
the slow rotators are indicated with open symbols. From 
Figure 1 we see that all three of the stars on which flares were 
seen have a dominant temperature component kT x > 1 keV. 
We also note that three of the four stars with k 7i < 1 keV are 
F stars, while the fourth is a slowly rotating G star. All stars 
with kT v > 1 keV are rapidly rotating G stars or K stars. No 
M stars had sufficient counts to constrain a two-temperature 
model. 

These results are in general agreement with those of Schmitt 
et al. (1990), who found generally hotter temperatures among 
K- and M-type dwarfs than among F- and G-type dwarfs 
observed by the Einstein IPC. We point out that their G star 
sample contained no stars with y sin i > 16 km s' \ That is, 
the G stars with '‘hot” coronae in the Pleiades are rotating 
much more rapidly than the G stars examined by Schmitt et al. 
(1990). Since most field G dwarfs are substantially older than 
their 70 Myr old counterparts in the Pleiades, the dearth of 
rapidly rotating field G dwarfs is not at all surprising given the 
timescales for stellar spin-down. 


In order to compare results for one- and two-temperature 
model fits, we introduce a mean temperature, 


<T X > 


L"«i ^em, 
ir=i EM,. ’ 


( 1 ) 


where N = 1,2 for one-temperature and two-temperature 
model results, respectively. 

For single-temperature models, 1 a upper and lower bounds 
were determined by fixing T and fitting EM for 0.1 
keV < kT <2.0 keV with 1 <x (68%) confidence corresponding 
1° X 2 ~ Xmin — 10. In order to determine confidence contours 
for the two-temperature models, a best-fit x 2 value was com- 
puted at each point in the kT { versus kT 2 plane by fitting the 
two remaining free parameters (EM t and EM 2 ). For two 
“interesting” parameters ( k 7\ and kT 2 ), 1 a (68%) confidence 
corresponds to x 2 — Xmin = 2.3 (Lampton, Margon, & Bowyer 
1976; Avni 1976). Using equation (1) and the 1 o minima and 
maxima of each temperature component, we have estimated 
the 1 a minima and maxima of </cT x >. 

In columns (1 1)— (13) of Table 1 we list <fcr x >, <kT min >, and 
</cT max ). Stars with no data in column (13) indicate that the 
upper bound is not constrained by the data. Also indicated in 
columns (14) and (15) are L x and L^JL^ L x has been com- 
puted by integrating the best-fit RS incident spectrum and 
assuming a cluster distance of 127 pc (Paper I) and a hydrogen 
column density of 2 x 10 20 cm" 2 . has been inferred from 
Table 7 of Paper I. 

We note that the X-ray luminosities have been determined 
using the best-fit incident spectrum and not by assuming a 
constant count-rate conversion factor, as was done in Paper L 
Consequently, the L x values in Table 1 are slightly different 
from those in Table 7 of Paper I. 

In Figure 2 we plot <fcT x > versus L^JL^ for the stars in 
Table 1. We use as a measure of X-ray activity in these 

late-type stars because this ratio allows us to estimate what 
fraction of the star’s total energy output is radiated in coronal 
soft X-rays. 

The most striking trend seen in Figure 2 and in Table 1 is a 
general increase in mean coronal temperature with X-ray 
activity. While the sample of stars for which X-ray spectra have 
been fitted is small and biased toward the most luminous stars, 



kT, (keV) 

Fig. 1.— Coronal temperatures of X-ray-bright late-type Pleiades stars 
based on RS two-temperature model fits. Plotted are kT x and kT 2 (keV) for F 
stars ( triangles ), G stars (circles), and K stars (squares), with open and filled 
symbols representing slow (i; sin i ^ 16 km s” ‘) and fast rotators, respectively. 
The error bars indicate l a uncertainties in kT t and kT 2 . 
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Fig. 2. — Mean coronal temperature vs. X-ray activity based on RS one- 
and two- tempera lure model fits. Plotted are </cT x > (keV) vs. log See 

Fig. 1 for symbol definitions. 
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this result, if correct, suggests that the coronae of the most 
active dwarfs do not simply contain more X-ray-emitting 
plasma than less active dwarfs, but that they contain generally 
hotter material as well. Again, a similar trend was noticed 
among late-type dwarfs, active binaries, and giants by Schmitt 
et al. (1990). 

As was demonstrated in Paper I, rotation is the dominant 
determinant of X-ray activity in the late-type Pleiades stars. 
Figure 1 suggests that lower effective temperature and higher 
rotational velocity may lead to higher coronal temperatures. 
From Figure 2, the trend is for increasing temperatures with 
increasing WUoi* If L \! L boi is a surrogate measure of the 
coronal heating efficiency, then we interpret Figure 2 as an 
increase in plasma temperature with increased heating. 
However, the small number of stars we have analyzed do not 
allow us to correlate coronal temperature directly with rota- 
tional velocity. 

Jordan & Montesinos (1991) have examined correlations 
between coronal emission measure, coronal temperature, and 
physical stellar parameters for a number of late-type dwarfs 
and giants detected by Einstein. Their results point to higher 
coronal temperatures with increasing activity and decreasing 
Rossby number. The coronal temperatures and volume emis- 
sion measures in Table 1, in fact, obey the same scaling laws. 
See Figure 3 of Jordan & Montesinos (1991). 

2.2. Composite PSPC Spectra of Cluster Stars 

While the results of the individual spectral fits may be sug- 
gestive, the few PSPC pulse-height spectra with sufficient 
counts to be well constrained are biased toward the most X- 
ray-luminous stars. In order to sample all the X-ray-emitting 
late-type Pleiades stars, we have generated composite spectra 
by co-adding spectra of individual stars according to spectral 
type and rotational velocity. Total source and background 
counts for all sources in each group were collected for each of 
the three exposures. The spectra were co-added by simply 
summing net counts and determining net rms errors in each 
energy channel. Since counts from individual sources have 


been extracted from different X-ray images and, hence, from 
different portions of the PSPC detector, the net energy 
response of the X-ray mirror assembly plus PSPC for the co- 
added spectrum has been calculated by determining a weighted 
mean of the individual energy responses. 

Since the composite spectrum is a simple sum of individual 
spectra of different sources, the composite spectrum may be 
dominated by counts from a handful of stronger sources. While 
this may bias results toward the stronger sources, an alterna- 
tive procedure (e.g., co-adding normalized spectra) would 
introduce larger net errors in the co-added spectrum, thereby 
poorly constraining the model parameters. 

The stars used in these composite spectra are the late-type 
cluster members [those with (B-F) 0 >0.35] in Table 7 of 
Paper I. In Table 2 we present a synopsis of the star samples 
used to generate the composite spectra. In column (1) we group 
the Pleiades stars according to spectral type and rotational 
velocity, and in column (2) we list the total number of stars in 
each group. As before, we define the rapid rotators as those 
stars with v vin i > 16 km s _1 . In column (3) we list the total 
net counts in the 0. 1-2.0 keV passband, and in columns (4)— (7) 
we list the best-fit kT l9 log (EM^EM*), and kT 2 , respec- 
tively, from the RS two-temperature model fits. 

As can be seen in column (3) of Table 2, for all but the slowly 
rotating F stars, > 1000 net counts have been collected in each 
composite spectrum, allowing relatively tight constraints to be 
placed on the model parameters. By the same procedure 
described in § 2.1, x 2 confidence contours were determined for 
the composite spectra in Table 2. 

In Figure 3 we present the x 2 68% and 90% confidence 
contours in the £7] versus kT 2 plane for all the F-type and 
K-type cluster members. Also plotted are the x 2 contours for 
rapidly and slowly rotating G stars. We have omitted the M 
stars for the sake of clarity. In actuality, the M-star contours 
are very close to those of the K stars. The crosses indicate the 
best-fit temperatures. While, the F-, K- and M-type composite 
spectra can be adequately fitted without regard to rotation, the 
spectrum made up of all G stars cannot be adequately fitted 


TABLE 2 

Composite ROSAT PSPC Spectral Fits 


Spectral Number Net xl Wi logEM 2 /EMi 

Group of Stars Counts (keV) (kcV) 

(1) (2) (3) (4) (5) (6) (7) 


F stars 

slow 4 601. 

fast 6 2668. 1.33 0.71 -0.26 0.15 

all 16 6381. 0.91 0.69 -0-22 0-17 


G stars 

slow 21 7890. 0.98 0.92 -0.33 0.23 

fast 12 14551. 1.03 1.27 -0.67 0.52 

all 33 22441. 1.83 — 


K stars 

slow.. . ... 24 5782. 1T26 L14 -0.28 0.30 

fast 16 7548. 1.14 1.08 -0.30 0.26 

all 47 13966. 1.35 1.11 -0.29 0.28 


M stars 

slow... ... 6 1233. 0^40 0X5 -0X4 0^20 

fast 5 1093. 0.68 1.15 -0.26 0.35 

aU 47 8678. 1.06 1.11 -0.31 0.26 
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Fig. 3. — Coronal temperatures of all X-ray-detected late-type Pleiades 
stars based on RS two- temperature model fits of the composite spectra. 
Plotted are the x 2 contours in the IcT, vs. kT 2 (keV) plane for F stars, slow G 
stars, fast G stars, and K. stars. The M star contours have not been indicated 
for the sake of clarity ; they are very close to the K. star contours. The inner and 
outer contour lines represent, respectively, the 1 a and 1.6 a (68% and 90%) 
confidence intervals in the fitted temperatures. 

with a two-temperature model. On the other hand, when we 
separate the fast and slow G stars, we find acceptable fits with 
substantially hotter distributions for the rapid rotators. We 
note that at the 90% confidence level, the model temperatures 
from the F, K, and rapidly and slowly rotating G stars do not 
overlap. 

Results of the RS model fits to the individual and composite 
spectra both point to increased coronal temperatures with 
decreasing effective temperature. There is also a suggestion 
that rotation may play a role in determining coronal tem- 
perature distributions, at least among the handful of ROSAT - 
detected G stars. This apparent dependence on rotation is 
probably an artifact of the overall increase in activity with 
rotational velocity. That is, the G-type Pleiades stars with 
higher rotational velocities have higher LyJL ^ ~ 10 3 and, 
consequently, hotter coronae. 

2.3. PSPC Hardness Ratios of Cluster Stars 

The coronal temperatures determined in the previous sec- 
tions are model-dependent, with different models (e g., one- 


temperature, two-temperature, or continuous emission 
measure distribution), different coronal abundances, and differ- 
ent plasma codes providing different best-fit plasma distribu- 
tions. An alternative approach to modeling the PSPC 
pulse-height spectra is to examine the hardness ratios of the 
detected sources. Analogous to optical photometric colours, 
X-ray hardness ratios have been used to estimate the level of 
circumstellar and interstellar hydrogen absorption. In addi- 
tion, accurate hardness ratio determinations can be used to 
estimate roughly the characteristic temperature of the emitting 
coronal plasma (e.g., see Fig. 4 of Neuhauser et al. 1995b). 

We use the standard ROSAT soft, medium, and hard bands 
defined, respectively, as s = net counts in channels 10-40 (0.1- 
0.28 keV), m = net counts in channels 51-90 (0.5-0.9 keV), and 
h = net counts in channels 91-200 (0.9-2.0 keV), and hardness 
ratios (Neuhauser et al. 1995a) as 


HR, 


h + m — s 


( 2 ) 


HR 2 


h — m 
h -f m 


( 3 ) 


We have determined 3 o count rates or upper limits in the 
three PSPC passbands for the 146 X-ray sources associated 
with a single F-, G-, K-, or M-type Pleiades cluster member. 
These count rates have yielded HR, and HR 2 values for 79 
stars. Because of upper limits in the soft band, in 58 cases only 
HR 2 could be determined. Because of upper limits in the 
medium or hard bands, in an additional nine cases neither 
HR, nor HR 2 could be determined. In Table 3 we present the 
hardness ratios and optical data for the 137 stars for which 
HR, or HR 2 could be determined. In columns (l)-(5) we list 
the X-ray source number, star name, (B— V) 0t (V -/) 0 , and v 
sin i from Table 7 of Paper L In columns (6)— (9) we list HR,, 
<t HR i , HR 2 and <7^. In columns (10) and (11) we list corre- 
sponding values of L x and Lx/L M from Paper I. 

In Figures 4 and 5 we plot the resulting hardness ratios 
versus LjJL^. See Figure 1 for symbol definitions. Figures 4 a 
and 4 b do not suggest higher HR, values with increased 
However, many of the K and M stars for which hard- 
ness ratios could be determined, i.e., the brightest ones, have 
values of W^boi close to the saturation value of 10 3 . The 
trend in Figure 5a is for increasing hardness ratio HR 2 with 



Fig. 4. — HR, vs. L x for (a) F stars ( trangles ) and G stars (circles) and (6) K 
i < 16 km s l ) and fast rotators, respectively. 



'°g LxAboi 

Fig. 4 h 

stars ( squares ) and M stars (diamonds). Open and filled symbols represent slow (v sin 




TABLE 3 

ROSAT PS PC Hardness Ratios 


X-ray 

Number 

(i) 


(B-V)o (V-I) o 


v suit 
(km s -1 ) 

( 5 ) 


2 

HCG 65 

4 

HCG 71 

6 

H II 26 

8 

H II 34 

9 

HCG 97 

10 

HCG 100 

12 

HCG 103 

17 

AKlB 121 

18 

H II 97 

24 

H II 134 

25 

HCG 115 

26 

H II 133 

30 

HCG 123 

31 

H II 173 

32 

H II 174 

35 

H II 191 

36 

H II 212 

37 

HCG 129 

39 

HCG 134 

40 

H II 253 

41 

H II 263 

42 

H II 250 

43 

SK 586 

46 

H II 293 

49 

H II 303 

53 

H II 314 

54 

H II 320 

55 

H II 324 

57 

H II 338 

58 

H II 335 

59 

HCG 149 

60 

HCG 148 

62 

H II 345 

63 

HCG 150 

64 

HCG 156 

65 

H II 347 

66 

H II 357 

69 

H II 405 

71 

H II 451 

72 

H II 476 

74 

H II 489 

79 

HCG 172 

80 

H II 514 

85 

H II 554 

89 

HCG 178 

90 

HCG 181 

92 

H II 625 

93 

H II 605 

94 

H II 624 

95 

H II 627 

97 

H II 686 

98 

H II 708 

100 

HCG 194 

102.... 

H II 738 
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H II 727 
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H II 745 

105.... 

H II 739 

107.... 

H II 746 


80 

0.63 

81 


31 

1.60 

33 

1.51 


1.96 

47 

2.03 

64 


84 

0.69 

.63 

0.45 


2.54 

.66 

0.49 

.84 

0.78 

.60 


.84 


.01 

1.03 

.42 


.22 

1.40 

.51 

1.97 




iO. 

47 

2.47 

* ■ • 


2.29 


37 

1.70 

75. 

17 

1.27 

<9. 

50 


17. 

17 

1.18 

<9. 

52 

0.22 

23. 

58 

0.44 

17. 


2.95 

. .. 

66 

0.48 

11. 

28 

1.42 


61 

2.50 

45 

53 

2.28 


.78 

0.55 

94. 

.40 


80. 

.49 

1.92 

16. 

.46 



.39 

1.32 

68. 

.58 


48 


2.37 


.75 

0.56 

50 

.51 


50 

.48 



.58 

0.44 

13 

.73 

0.61 

<7 


-3.20 

-3.21 

-4.89 

-4.04 

- 2.88 

-2.87 

-3.26 

-3.33 

-3.46 

-3.16 

- 2.68 

-3.01 

-2.96 

-4.25 

-2.93 

-2.98 

-3.04 

-3.11 

- 2.88 

-3.01 

-3.33 

-4.06 

-2.83 

-4.29 

-3.72 

-3.22 

-3.47 

-2.83 

-4.77 

-2.89 

-2.82 

-2.58 

-2.98 

-2.65 

-3.28 

-3.12 

-2.83 

-4.13 

-3.42 

-4.17 

-4.24 

-2.62 

-4.22 

-3.31 

-2.77 

-2.52 

-2.92 

-5.16 

-3.20 

-4.34 

-2.93 

-3.89 

-2.71 

-3.04 

-3.85 

-4.41 

-3.61 

-4.17 
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(B - V) 0 

0 

1 

v sin i 

HR\ 


hr 7 

<r 2 

Lx 

^bol 

Number 

Name 



(km a -1 ) 





(erg» 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

108 

H II 762 

1.36 

1.48 

<9. 

. . . 


0.01 

0.10 

29.33 

-3.16 

109 

H II 761 

0.63 


11. 

0.33 

0.05 

-0.11 

0.05 

29.85 

-3.66 

116 

SK 488 


2.10 


0.64 

0.12 

0.08 

0.08 

29.25 

-2.85 

116 

H II 885 

0.98 

0.96 

<8. 

... 


0.08 

0.26 

28.98 

-4.05 

117 

H II 890 

1.29 

1.66 

<9. 

0.38 

0.18 

0.24 

0.18 

29.02 

-3.34 

119 

H II 906 

1.48 

2.19 


0.49 

0.13 

-0.03 

0.08 

29.51 

-2.87 

120 

H n 930 

1.22 

1.48 

20. 

0.52 

0.14 

-0.06 

0.12 

29.15 

-3.20 

121 

H II 916 

0.83 


<7. 

0.07 

0.12 

0.08 

0.12 

29.22 

-3.87 

126 

H II 996 

0.59 

0.41 

10. 

0.55 

0.08 

0.03 

0.07 

29.55 

-4.01 

130 

HCG 219 

1.55 

2.26 


0.47 

0.17 

0.06 

0.15 

29.06 

-3.09 

132 

H II 1039 

0.93 

0.66 

<7. 

... 

... 

0.08 

0.11 

29.85 

-3.14 

134 

H II 1032 

0.71 

0.57 

37. 

0.73 

0.03 

0.10 

0.03 

30.16 

-3.14 

136 

H II 1061 

1.15 

1.23 

<9. 

0.72 

0.12 

0.14 

0.09 

29.51 

-3.01 

141 

H II 1094 

1.36 

1.56 



... 

0.01 

0.09 

29.60 

-3.00 

142 

H II 1100 

1.06 

1.02 

<9. 

0.66 

0.05 

0.07 

0.04 

29.87 

-3.14 

143 

H II 1117 

0.69 


<7. 

0.26 

0.12 

-0.02 

0.13 

29.33 

-4.33 

144 

H II 1122 

0.42 


... 

0.59 

0.07 

-0.18 

0.06 

29.55 

-4.43 

145 

H II 1136 

0.72 

0.57 

75. 


... 

0.08 

0.07 

30.14 

-3.02 

146 

H II 1124 

0.94 

0.81 

<8. 

0.56 

0.10 

0.06 

0.08 

29.47 

-3.43 

147 

H II 1101 

0.57 

0.42 

18. 


... 

-0.27 

0.18 

29.37 

-4.24 

163 

H II 1215 

0.58 

0.42 

<7. 


... 

-0.20 

0.16 

29.34 

-4.16 

154 

HCG 244 


2.45 

... 


... 

-0.05 

0.11 

29.21 

-2.91 

156 

HCG 246 


2.67 

... 

0.22 

0.15 

-0.01 

0.17 

29.11 

-2.60 

160 

H II 1280 

1.40 

1.61 

75. 

0.62 

0.08 

0.05 

0.07 

29.44 

-2.65 

161 

H II 1298 

0.96 

0.82 

<7. 

0.66 

0.11 

-0.04 

0.08 

29.51 

-3.40 

164 

H II 1321 

1.43 

2.17 

17. 


... 

-0.37 

0.19 

29.16 

-3.20 

167 

H II 1309 

0.43 



0.51 

0.06 

-0.21 

0.05 

29.60 

-4.32 

168 

HCG 258 

1.50 

2.27 




-0.01 

0.20 

29.27 

-2.91 

169 

H II 1332 

0.98 

0.84 

<8. 

0.37 

0.08 

-0.02 

0.09 

29.57 

-3.28 

170 

H II 1338 

0.42 


10. 

0.42 

0.13 

0.00 

0.12 

29.12 

-5.10 

172 

H II 1365 

1.36 

1.65 

<9. 

0.42 

0.07 

-0.03 

0.07 

29.59 

-3.04 

174 

H II 1348 

1.11 

1.12 

<9. 


... 

-0.08 

0.16 

29.09 

-3.77 

181 

H II 1454 

1.04 

0.98 

<9. 


... 

-0.20 

0.21 

29.16 

-3.59 

182 

HCG 273 


1.86 


0.42 

0.13 

0.06 

0.10 

29.26 

-2.83 

187 

H II 1516 

1.30 

1.43 

25. 

0.60 

0.05 

0.11 

0.05 

29.89 

-2.64 

188 

H II 1514 

0.62 

0.44 

15. 

0.61 

0.09 

0.03 

0.07 

29.27 

-4.24 

189 

H II 1632 

1.23 

1.37 

32. 

0.38 

0.13 

-0.03 

0.12 

29.43 

-3.01 

190 

H II 1531 

1.23 

1.22 

50. 

0.45 

0.08 

0.00 

0.08 

29.56 

-3.04 

195 

H II 1613 

0.50 


19. 

0.24 

0.12 

-0.59 

0.15 

29.30 

-4.46 

197 

H II 1653 

1.17 

1.25 

35. 

0.31 

0.15 

0.00 

0.09 

29.59 

-2.99 

198 

H II 1726 

0.51 


12. 

0.61 

0.08 

-0.21 

0.06 

29.56 

-4.45 

210 

H II 1827 

1.46 

1.96 


... 


-0.01 

0.18 

29.28 

-3.15 

211 

H n 1856 

0.52 


16. 

... 


-0.15 

0.11 

29.10 

-4.60 

217 

H II 1912 

0.45 


75. 

0.58 

0.10 

0.08 

0.07 

29.43 

-4.65 

219 

HCG 337 


2.37 


... 

... 

-0.07 

0.25 

29.05 

-2.81 

222 

H II 2034 

0.92 

0.89 

77. 

0.59 

0.08 

-0.05 

0.07 

29.67 

-3.09 

223 

H II 2027 

0.81 

0,65 

<8. 

0.45 

0.12 

0.03 

0.09 

29.42 

-3.98 

224 

HCG 349 


2.48 


0.34 

0.15 

0.03 

0.11 

29.26 

-2.93 

228 

H II 2147 

0.76 


27. 

0.62 

0.02 

0.11 

0.02 

30.60 

-2.82 

229 

HCG 355 


2.70 


0.14 

0.14 

0.03 

0.13 

28.83 

-2.85 

231 

H II 2172 

0.58 

0.39 

10. 


... 

-0.21 

0.17 

29.30 

-4.26 

234 

H II 2208 

1.32 

1.64 

73. 

0.27 

0.12 

0.21 

0.09 

29.42 

-3.05 

237 

H II 2244 

0.95 

0.92 

45. 

0.76 

0.07 

0.16 

0.05 

29.99 

-2.82 

241 

H II 2278 

0.83 

0.67 

<9. 

... 


0.20 

0.17 

29.58 

-3.82 

243 

HCG 370 


2.24 


... 

... 

0.04 

0.20 

28.82 

-3.29 

247 

HCG 372 

1.48 

2.41 


0.40 

0.17 

-0.09 

0.17 

29.05 

-2.87 

249 

HCG 376 


2.48 


0.33 

0.11 

0.00 

0.12 

29.00 

-2.69 

260 

H II 2366 

0.78 

0.59 

<7. 

0.33 

0.11 

0.00 

0.12 

28.85 

-4.32 
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X-ray 

Star 

(B - V) 0 

(v - /) o 

v sint 

HRi 

0\ 

hr 2 

<7 2 


L x 

Number 

Name 



(km s' 1 ) 





(ergs s" 1 ) 


(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(8) 

(10) 

(H) 

253 

H II 2407 

0.91 

0.75 

<7. 


... 

0.07 

0.14 

29.16 

-3.78 

266 

H II 2548 

1.29 

1.39 

<9. 



0.04 

0.20 

28.98 

-3.52 

268 

H II 2602 

1.59 

2.25 

11. 

0.53 

0.16 

0.04 

0.15 

29.42 

-2.87 

269 

H II 2601 

1.49 

2.12 

10. 

0.42 

0.11 

0.03 

0.11 

29.37 

-3.05 

270 

H II 2588 

1.14 

1.13 

<7. 

0.44 

0.16 

0.21 

0.15 

29.20 

-3.42 

272 

HCG 394 

1.93 

2.74 

... 

... 


0.07 

0.23 

28.95 

-2.74 

274 

H II 2644 

0.73 

0.49 

<7. 

... 


-0.06 

0.18 

29.04 

-4.27 

280 

H II 2741 

0.97 

0.87 

10. 



0.02 

0.10 

29.38 

-3.41 

282 

H II 2786 

0.56 


25. 

... 


-0.12 

0.15 

29.38 

-4.21 

287 

H II 2881 

0.93 

0.84 

12. 

0.50 

0.09 

0.01 

0.08 

29.90 

-3.30 

289 

H II 2880 

0.82 


<7. 


... 

-0.20 

0.24 

28.98 

-4.09 

292 

H II 2927 

1.22 

1.37 

95. 



-0.08 

0.16 

29.34 

-3.10 

293 

HCG 422 


2.37 


0.55 

0.14 

0.17 

0.12 

29.30 

-2.78 

294 

HCG 424 


2.31 


0.34 

0.18 

-0.06 

0.18 

29.22 

-2.85 

295 

H II 3019 

1.17 

1.20 

<9. 



-0.03 

0.17 

29.28 

-3.29 

296 

HCG 428 


2.04 

... 

... 


-0.11 

0.16 

29.43 

-2.63 

297 

H II 3063 

1.13 

1.23 

28. 

0.42 

0.16 

0.09 

0.14 

29.61 

-2.93 

298 

H II 3096 

0.89 

0.76 


0.51 

0.08 

-0.06 

0.07 

29.74 

-3.23 

300 

H II 3163 

0.97 

0.97 

70. 

0.50 

0.08 

-0.02 

0.07 

29.78 

-2.96 

302 

H II 3179 

0.52 


<7. 


... 

-0.17 

0.18 

29.27 

-4.42 

304 

H II 3197 

1.03 

1.07 

33. 

0.54 

0.05 

0.03 

0.05 

30.14 

-2.91 


increasing Lx/L^, for the F- and G-type stars. More signifi- 
cantly, the rapidly rotating G stars all like in a quadrant 
defined by HR 2 > 0.0 and > -3.3, while the F stars 

and slowly rotating G stars generally lie in the opposite quad- 
rant. If we interpret LyJL^ as the efficiency of the coronal 
heating mechanism, then Figure 5a suggests increased coronal 
temperatures with increased coronal heating efficiency. On the 
other hand, there is little we can say about the K- and M-type 
hardness ratios because the relative errors are so large and 
because we are only sampling the most active lower mass stars. 

A formal analysis of the error-weighted correlations (see 
Bevington & Robinson 1992) indicates that HR 2 is correlated 
with L x and W^i at the 99% confidence level for the F- and 
G-type Pleiades stars. For the K and M stars, HR t is corre- 
lated with L x at the 99% confidence level. These stars have 



hardness ratios and activity levels more comparable to those of 
the rapidly rotating G stars. 

In conclusion, the analysis of the individual X-ray spectra, 
the composite X-ray spectra, and the hardness ratios derived 
for the late-type stars all point to cooler coronal temperatures 
and lower activity levels for F stars and slowly rotating G-type 
dwarfs when compared to those of the K-, M-, and more 
rapidly rotating G-type dwarfs. 

As has been noted, these results are similar to those found by 
Schmitt et al. (1990) for late-type field stars. When comparing 
overall X-ray activity from field stars and open clusters, L x has 
been shown, overall, to decrease with age (e.g., Rosner, Golub, 
& Vaiana 1985). The similarity, then, between the coronal dis- 
tributions in the Pleiades and the stars examined by Schmitt et 
al. (1990) may appear somewhat surprising. We point out, 


0.4b 



-0.4b 


- 0.61 ‘ - 

-5.0 -4.5 -4.0 -3.5 -3.0 -2.5 

log L x /L bol 

Fig. 5b 

K stare and M stars. See Fig. 4 for symbol definitions. 
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however, that the Schmitt et al. sample is biased toward 
extremely bright X-ray sources, which, on average, have X-ray 
activity levels more comparable to those of young dwarfs like 
the Pleiades than of typically less active field dwarfs. 

3. X-RAY FLARES 

The ROSAT PSPC data contain the time of arrival of every 
detected X-ray photon. From the photon event data we have 
constructed background-subtracted light curves of every 
detected source in every one of the long exposures in Table 1 of 
Paper I, allowing us to look for source variability on timescales 
as short as 500 s. 

Each light curve was visually inspected for signs of variabil- 
ity or obvious fiaring. We also used a simple x 2 test on each 
light curve to test the null hypothesis of a constant count rate. 
While a number of sources were found not to be constant 
(above the 99% confidence level), only half or so of these 
showed the characteristic rise and decay of a strong stellar 
flare. We found no sources showing strong evidence of periodic 
variability. 

In Figure 6 we present light curves for the 12 most convinc- 
ing cases of X-ray fiaring in the PSPC data. In Figure 6 we plot 
the observed count rate (not corrected for telescope vignetting) 
with 1 <7 error bars against time from the beginning of the 
observation. Ten of these flares are associated with known, 
late-type Pleiades cluster members. The most intriguing of the 
flares is associated with source 191. This source has no cata- 
loged stellar counterpart; however, a very faint star (K>16) 
appears less than 10" away on the Palomar Observatory Sky 
Survey E plate of the region (see Fig. 8 of Paper I). If the 
counterpart to the observed flare is indeed a Pleiades cluster 
member, then the peak X-ray luminosity during flare would 
have reached ~ 1 0 3 1 ergs s “ 1 . 

All of the flares associated with Pleiades cluster members 
had peak flare luminosities in excess of 10 30 ergs s' 1 , the 
brightest one, on the G star H II 2147, with L£* k ~~ 2.5 x 10 31 
ergs s “ 1 . The flares are generally characterized by rise times in 
the range (1-10) x 10 3 s, decay times in the range (5-15) x 10 3 


s, and decay times typically exceeding rise times by a factor of 
-3. 

In Table 4 we summarize our findings for the dozen flares in 
Figure 6. In columns (1)— (4) we list the X-ray number, star 
name, ( B — K) 0 , and (V — I) 0 colors of the optical counterpart. 
In columns (5)-(8) we list the quiescent and peak X-ray lumi- 
nosities as well as the e-folding times of the rise and decay 
phases of the flare. Only two flare spectra contained sufficient 
counts to model the peak of the emission: H II 2147 and H II 
1516 (see § 3.1 below). 

For the sake of consistency, we have assumed the same 
count rate to luminosity conversion factor as was used in 
Paper I for all the flares in Table 4. Since stellar flares are 
nearly always characterized by temperatures in excess of 1 keV 
(i.e., >12 x 10 6 K) (Poletto, Pallavicini, & Kopp 1988; van 
den Oord & Mewe 1989; Fisher & Hawley 1990; Schmitt 
1994), we have probably underestimated the true count-rate to 
luminosity conversion factor. Consequently, the L%* k values in 
Table 4 are lower limits. The rise-time and decay-time upper 
limits in Table 4 indicate that the rise or the decay phase of the 
flare was not sufficiently monitored to determine an e-folding 
time. 

Overall, the Pleiades flares are among the brightest X-ray 
flares ever reported on single, late-type main-sequence dwarfs. 
Stellar X-ray flares with L x ~ 10 31 ergs s _l have only been 
reported for the eclipsing binary /? Per (Algol) (van den Oord & 
Mewe 1989), for the flare star AU Mic (Cully et al. 1994), and 
for a handful of pre-main-sequence stars in p Oph (Montmerle 
et al. 1983) and Orion (Gagne et al. 1995) and for the Pleiades 
star H II 1136 as observed by the Einstein IPC (Caillault & 
Helfand 1985). During the ROSAT PSPC pointed observation, 
H II 1136 was not variable above the 99% confidence limit, 
and no flare was apparent. 

Comprehensive surveys of Einstein and EX OS AT obser- 
vations of stellar X-ray flares have been presented by Haisch 
(1983) and Pallavicini, Tagliaferri, & Stella (1990), respectively. 
While most of those X-ray flares occurred on dMe stars, the 
ROSAT all-sky survey has found that large flares occur on all 
types of late-type stars (Schmitt 1994). Our results confirm the 
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Number 

Name 



Quiet 

Peak 

Rise 

Decay 

(i) 

(2) 

(3) 

(4) 

(5) 

(«) 

(7) 

(«) 

9 

HCG 97 


2.02 

29.2 

30.2 

<3.4 

5.0 

32 

H II 174 

0.81 

... 

30.0 

30.5 

5.5 

13.6 

35 

H II 191 

1.31 

1.61 

29.3 

30.5 

1.0 

5.0 

36 

H II 212 

1.33 

1.51 

29.3 

30.3 

<4.0 

6.0 

56 

HCG 143 

1.61 

2.37 

29.6 

30.4 4 

0.5 

1.5 

62 

H II 345 

0.80 

... 

30.1 

30.7 

<4.0 

6.5 

90 

HCG 181 

1.53 

2.28 

29.5 

30.5 

1.0 

5.0 

142 

H II 1100 

1.06 

1.02 

29.7 

30.5 b 

3.0 

9.5 

187 

H II 1516 

1.30 

1.43 

29.9 

31.2 

<5.0 

5.0 

191 




29.7 

30. 9 C 

<4.5 

11.0 

228 

H II 2147 

0.76 


30.6 

31.4 b 

<11.0 

<84.0 

237 

H II 2244 

0.95 

0.92 

29.9 

30.6 

<5.0 

8.0 


* Although the membership of HCG 143 is uncertain, we assume a distance of 127 pc. 
b The peak of the flare was probably not observed. 

c Although this source has no optical counterpart, we assume a distance of 127 pc. 
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Fig. 6e 

Fig. 6.— Light curves for the dozen strong X-ray flares detected in the ROSAT PSPC data of the Pleiades. Plotted are the uncorrected PSPC count rates vs. time 
in seconds from the start of the observation. The error bars indicate 1 a uncertainties in the count rate. The light curves have bin lengths of 1000 s. Also indicated are 
the field from which the light curve was extracted, the observation start date, and the X-ray source number (see Table 4). 


all-sky survey finding with flares detected on dwarfs ranging in 
spectral type from -G8 (H II 2147) to - M3 (HCG 143) and 
ranging in rotational velocity from t; sin i < 9 km s 1 (H II 
1 100) to —45 km s~ 1 (H II 2244). 

Although our sample of detected X-ray flares is small and 
strongly biased toward very powerful events, we can begin to 
estimate the flaring frequency of the Pleiades dwarfs. The large 
flares we detect had typical durations (rise + decay) of 10 4 s 


with typical monitoring periods of — 6 x 10 4 s. Of the 171 
detected cluster members, 10 were seen to flare. We conclude 
that the flaring frequency for large (L x > 10 30 ergs s' J ), long- 
duration (t - 10 4 s) flares is about 1%. This appears to be 
higher than the flaring frequency observed in the Hyades 
(Stern et al. 1994) and may be indicative of the comparatively 
lower activity levels of the Hyades dwarfs. This decreased 
X-ray activity is presumably the result of the overall decrease 
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in rotation rates from the age of the Pleiades (t* — 70 Myr) to 
the age of the Hyades (r + ~ 600 Myr). 

3.1. Modeling the Flares on H II 2147 and H II 1516 

Two of the detected flares had sufficiently high count rates to 
model the PSPC spectrum of the flare peak. In the case of H II 
2147, the flare was only seen during one ROSAT orbit with a 
total exposure time of ~ 1900 s. The ^ 1000 net counts in the 
flare spectrum were modeled with a two-temperature RS 
model, the results of which we present in Table 1. Since the rise 
and decay of the flare were not observed, there is little else we 


can infer from the H II 2147 light curve. We point out the 
enormous emission measure, EM > 2 x 10 54 cm' 3 , during 
flare peak. 

The second brightest flare is associated with the late-K type 
Pleiades star H II 1516. We estimate the e-folding time of the 
decay phase to be ^5.0 x 10 3 s. The peak of the flare was 
observed for ~660 s during one ROSAT orbit. The ^400 net 
counts during this orbit were modeled with a one-temperature 
RS model. Although the signal-to-noise ratio of the PSPC 
spectrum of the flare peak is very low, the minimum flare 
temperature and emission measure we obtain (at the 68% con- 
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fidence level) are kT >1.1 keV and 1.2 x 10 54 cm 3 < EM < 
2.4 x 10 54 cm -3 . The flare temperature (T £ 13 MK) is typical 
of large flares on late-type dwarfs. We cannot constrain the 
temperature of a very hot plasma because of the low signal-to- 
noise ratio of the X-ray spectrum and the low sensitivity of the 
PSPC at energies above ~ 1.4 keV. As a result, we estimate the 
peak temperature to be T 0 ~ 25 MK, a typical value for late- 
type flares (Schmitt 1994; Pallavicini et al. 1990). As was the 
case for H II 2147, the volume emission measure of H II 1516 
at the onset of the flare, EM ~ 2 x 10 54 cm “ 3 , is very large. 

We can use the observed e-folding time of the decay phase to 
estimate the average electron density and the total volume of 
the plasma at the onset of the flare. In the quasi-static cooling 
loop model (van den Oord & Mewe 1989), one assumes that 
the flaring plasma decays through a series of stationary loops 
at constant volume and pressure. These assumptions are in 
contrast to the modelling schemes of Poletto et al. (1988) and 
Fisher & Hawley (1990). However, since the latter two models 
are more complex and, hence, require more detailed knowledge 
of the flare decay, we will limit our analysis to quasi-static 
cooling. 

We follow Schmitt (1994) and define two characteristic time- 
scales for the decay phase : the observed e-folding time, ij, and 
the radiative cooling time defined as 

(4) 

Tc 

where T 0 , n 0 , and A(T 0 ) are the initial plasma temperature, 
mean electron density, and radiative cooling law, respectively, 
at the beginning of the decay phase. In the quasi-static loop 
model, if one assumes no additional heating during the decay 
phase, the time evolution of the total radiative output reduces 
to a particularly simple form : 

E r Jt) = 4 + • (5 > 

At t = t d , E,Jr d )!E 0 = e~\ Consequently, we can infer the 
radiative cooling time directly from the decay time, 




Tc 3(e(l/4) 1/2 - 1) 


1.2 T d . 


( 6 ) 


In the case of H II 1516, x d - 5.0 x 10 3 s. Using an estimated 
temperature T 0 ~ 2.5 x 10 7 K, assuming a cooling law appro- 
priate for this temperature regime A(T 0 ) ~ 1.3 x 10 23 ergs 
cm 3 s _1 (van den Oord & Mewe 1989), substituting for x d in 
equation (6), and solving for equation (4), we obtain a mean 
electron density n 0 ~ 1.3 x 10 11 cm 3 . This is certainly a 
lower limit to the peak electron density because the true peak 
plasma temperature may, in fact, be higher than our estimate. 

If we assume a constant electron density throughout the 
flare loop, we infer from the peak observed emission measure, 
EM — 2 x 10 54 cm 3 , a total plasma volume V ~ 10 32 cm 3 . If 
the flare is the result of a single coronal loop with aspect ratio 
(ratio of loop diameter to loop length) a - 0. 1 (see van den 
Oord & Mewe 1989), then the inferred loop size is huge, 
L - 3 x 10 1 1 cm. That is, L > 6 R*. Of course, the flare may be 
the result of multiple loops, in which case the loop length can 
be scaled down accordingly. 

We point out, though, that the large plasma volume may be 
the result of our underestimate of the peak flare temperature or 
may be an artifact of the quasi-static cooling loop model. As 


has been pointed out by Schmitt (1994), alternative models like 
the two-ribbon flare model of Poletto et al. (1988) lead to 
higher electron densities and more compact loop structures. 
Higher energy X-ray spectra and more complete temporal 
coverage of flare events are required to model X-ray flares 
correctly. 

4. MEDIUM- AND LONG-TERM X-RAY VARIABILITY 

The PSPC data for the Pleiades were collected in exposures 
separated by 6-18 months. The Pleiades was also monitored 
by Einstein more than a decade ago, allowing us to look at a 
particular source’s variability on medium- and long-term time- 
scales. 

Throughout this section we employ a procedure for estimat- 
ing the cumulative count-rate ratio distribution given a pair of 
observations of a set of X-ray sources first used by Schmitt et 
al. (1993) and subsequently by Gagne & Caillault (1994) and 
Gagne et al. (1995). 

For each pair of observations of a particular source, we 
compute a count-rate ratio, r = w h ere an ^ are the 

count rates or upper limits from the earlier and later obser- 
vation, respectively. Since some sources detected in one expo- 
sure are observed but not detected in others, the count-rate 
ratios consist of some upper and lower limits. 

We bin the three sets of log r values into three distributions 
(count-rate ratios, upper limits, and lower limits). The set of 
r-values for all observations are doubly censored univariate 
data, and we estimate the “true” distribution of count-rate 
ratios iteratively using the maximum likelihood estimator 
described by Schmitt et al. (1993) 

4.1. Medium-Term Variability: 1 Year Timescale 

The PSPC exposures of the northeast and northwest fields 
were all observed in 1991 September. However, the obser- 
vations of the center field were done in three segments: 1991 
February, 1991 August, and 1992 August with exposure times 
of ^16, and ks, respectively. This has allowed us to 
look at count-rate variations on three timescales: 6 months 
between 1991 February and 1991 August, 12 months between 
1991 August and 1992 August, and 18 months between 1991 
February and 1992 August. To eliminate any bias due to the 
different exposure times of the three observations, we have 
further divided the individual exposures into appropriately 
long blocks. That is, to compare the 1991 February and 1991 
August data sets, the 1991 August data have been divided into 
four contiguous blocks of ^ 4 ks. 

For each X-ray source associated with a single Pleiades 
cluster member, we have determined a count-rate or 3 cr upper 
limit in each block of time. For each pair of observations, a 
count-rate ratio has been determined. Because of the presence 
of upper limits in both the “before” and “after data, the 
count-rate ratios consist of ratios, upper limits, and lower 
limits. Ratios determined from two upper limits have been 
discarded. We have used the iterative maximum likelihood 
technique described by Schmitt et al. (1993) and Gagne & 
Caillault (1994) to estimate the true count-rate ratio distribu- 
tion. 

In all cases, the distributions are symmetric about log r - 0 
and, hence, have been folded about log r = 0 to construct 
cumulative distributions. In Figure 7 we present the resulting 
cumulative count-rate ratio distributions from observations 
separated by 18 months ( top panel), 12 months ( middle panel), 
and 6 months (bottom panel) as solid histograms. 
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Fig. 7.—Count-rate ratio histograms for sources associated with late-type 
Pleiades stars from observations separated by 18 months ( top panels 12 
months (middle panef), and 6 months {bottom panel). The solid histograms 
represent the cumulative maximum likelihood count-rate ratio distribution for 
detections and upper limits using the estimator for doubly biased univariate 
data described by Schmitt et al. (1993). The dotted lines indicate the cumula- 
tive count-rate distribution we expect as the result of statistical fluctuations 
alone in the observed data in the case of a faint source with an intrinsic mean 
count rate at the 3 o detection threshold. 

If the X-ray count rates were very precise, then the resulting 
distribution of count-rate ratios would indicate intrinsic X-ray 
variability. However, since many of the sources are detected at 
or near the source detection threshold, repeated measurements 
of an intrinsically constant source will produce a broad dis- 
tribution of count-rate ratios as the result of measurement 
errors alone. We want to estimate how much of the observed 
variability is intrinsic and how much is the result of statistical 
fluctuations alone. Although the fluctuations from a photon 
counter like the PSPC are, strictly speaking, Poisson distrib- 
uted, all the sources in this survey have been detected with at 
least 20 net photons; thus, we assume that a Poisson error 
distribution can be accurately approximated by an appropriate 
Gaussian distribution. 

We consider the case of a source with a constant intrinsic 
count rate near the 3 a source detection threshold, where a is 
the total uncertainty in the source and background counts (see 
Gagne & Caillault 1994). In this case, the probability of a 
measurement yielding an upper limit is 0.5, and, after repeated 
measurements, variations due to statistical fluctuations alone 
will result in count-rate ratios as well as some upper and lower 
limits. The resulting distribution of count-rate ratios, upper 


limits, and lower limits is used to calculate a maximum likeli- 
hood distribution as before. 

In Figures 7, 8, and 9 (see below) we indicate the estimated 
cumulative distribution due to “errors only ” with dotted lines. 
One of the important consequences of such an error analysis is 
that, because of the low signal-to-noise ratio of the detected 
sources, count-rate variations of less than a factor of 2 may not 
be the result of intrinsic variability but, rather, may be the 
result of measurement uncertainties. 

In Figure 7 the difference between the solid histograms and 
dotted curves represents a lower limit to the true variability of 
the X-ray sources. As an ensemble, £ £ of the Pleiades stars are 
variable by more than a factor of 2 on all three timescales. 
Since the three distributions are so similar, we have combined 
the data to increase our sample size. 

In Figure 8 we present the combined cumulative count-rate 
ratio histogram indicating the observed variability on 6-18 
month timescales. The error bar is a typical 1 a uncertainty in 
each log r bin. Uncertainties in the cumulative distribution 
were determined by adding and subtracting 1 a uncertainties in 
the individual count rates, redetermining a set of count-rate 
ratios (with upper and lower limits) and recomputing the 
cumulative maximum likelihood distribution in each case. 

Since only two strong flares (H II 2147 and H II 1100) were 
identified in the center field, the variability observed in over 
one-quarter of the observations is not only the result of very 
strong X-ray flaring. Since we have not identified more moder- 
ate flaring (L x ^ 10 3 ° ergs s " *), we can only guess at the origin 
of this variability. We can eliminate fluctuations in the PSPC 
detector sensitivity as a possible explanation because such 
variations would have affected all the count rates determined 
at a particular time. This, in turn, would have skewed the log r 
values, giving rise to an asymmetric count-rate ratio distribu- 
tion. This was not observed. We can also eliminate any uncer- 
tainties in the off-axis response of the PSPC because all the 
count rates were determined from the center field alone. That 
is, count rates for a given source were all collected from the 
same position on the detector using identical source and back- 
ground cells. 

Our finding that about one-fourth of the late-type Pleiades 
stars vary by more than a factor of 2 on timescales ~ 1 yr is in 



Fig. 8. — Combined cumulative count-rate ratio histogram for sources 
associated with late-type cluster members. Overall, one-fourth or more of the 
Pleiades stars are significantly variable (by more than a factor of 2) on 6-18 
month timescales 
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contrast to the results of Micela et al. (1990, 1995), who found 
little significant variability in their analyses of the entire Ein- 
stein data set and a single ROSAT PSPC pointing of the 
southwest portion of the cluster. 

On the other hand, the fraction of variable Pleiades dwarfs 
presented herein is very similar to that observed among the 
very active pre-main-sequence stars in the Orion Nebula 
region by Gagne & Caillault (1994) and Gagne et al. (1995). In 
a qualitative sense, the variability in the Pleiades is also remi- 
niscent of the X-ray luminosity variations seen in the p Oph 
star-forming region (Montmerle et al. 1983), albeit on a smaller 
scale. 

4.2. Long-Term Variability: 11 Year Timescale 

The Pleiades open cluster was the object of intense monitor- 
ing by the Einstein IPC, the most comprehensive results of 
which have been reported by Micela et al. (1990). In their 
analysis of the ROSAT all-sky survey data of the Pleiades, 
Schmitt et al. (1993) found that ~40% of the sources were 
variable by more than a factor of 2 in the 11 yr time span 
between the all-sky and Einstein observations, and suggested 
that the variability and the decade-long time span between 
observations may be indicative of solar-like activity cycles. 
Solar-like activity cycles have been observed in some late-type 
dwarfs from monitoring of chromspheric Ca ii and Mg n flux 
variations and from optical continuum variations due to spot 
activity (Wilson 1978; Dorren & Guinan 1994; Baliunas et al. 
1995). 

Although the Einstein observations of the Pleiades had 
varying exposure times, X-ray luminosities were typically 
L x > 1.5 x 10 29 ergs s' 1 , a factor of 2 or 3 less sensitive than 
the ROSAT pointed observations. In order to avoid possibly 
biasing our analysis, we consider only cluster members that 
could have been detected in both sets of observations Of the 
detected sources reported by Micela et al. (1990), 69 were 
observed in one of the three pointed ROSAT observations. 

We thus consider the 69 brightest Pleiades stars detected by 
ROSAT and observed by Einstein . Of these, 51 were detected 
by Einstein and 18 had upper limits. Of the remaining 18 Ein- 
stein detections, 14 were detected by ROSAT and four yielded 
only ROSAT upper limits (see Table 4 of Micela et al. 1990 and 
Tables 7 and 8 of Paper I). We thus have 87 
(51 4- 18 4 14 4 4 = 87) pairs of observations consisting of 65 
Pleiades cluster members detected in both surveys, four 
members detected only by Einstein , and 18 members detected 
only by ROSAT. The resulting X-ray luminosity ratios consist 
of ratios, upper limits, and lower limits. 

Although the PSPC and the IPC have slightly different 
bandpasses, the median value of log L%** ei */LZ° SAT is approx- 
imately 0.0, indicating overall consistency in the soft X-ray 
luminosity determinations in the two surveys. We have calcu- 
lated the cumulative, maximum likelihood distribution as 
before. 

In Figure 9 we present the resulting cumulative L x ratio 
distribution ( solid line) and the distribution we estimate due to 
measurement uncertainties alone (dotted line). The error bar is 
a typical 1 o uncertainty in the ratio distribution itself. 
Whereas ~ 25% of the Pleiades stars are variable by more than 
a factor of 2 on 6-18 month timescales, up to 40% are variable 
on the — 11 yr timescale between Einstein and ROSAT. Similar 
variability was observed by Schmitt et al. (1993) in their 
analysis of the ROSA T all-sky data. 

We note that the sample used to compare the Einstein and 



Fig. 9. — Count-rate histogram for the 69 brightest sources observed by 
both the K0S/4T PSPC and the Einstein IPC. Approximately 40% of the 
Pleiades stars are significantly variable (by more than a factor of 2) on the ~ 1 1 
yr timescale between missions. 

ROSAT pointed observations are brighter in X-rays than the 
sample used to examine approximately 1 yr variability in the 
ROSAT data alone. We have examined the effect of this bias by 
recomputing count-rate distributions (as was done in § 4.1) for 
a reduced sample of ROSA T sources with the same flux limit 
as the Einstein observations, i.e., corresponding to L x > 1.5 
x 10 29 ergs s “ l . We find, as before, that approximately 25% of 
sources are variable by more than a factor of 2 on roughly 1 yr 
timescales. We conclude that we are seeing marginal evidence 
for increased activity over the decade-long timescale. 

The fraction of X-ray sources in the Pleiades that do show 
some variability (around 25%) do so on all timescales longer 
L han a month. Hence, any solar-like activity-cycle signatures 
would manifest themselves as additional variability. The 
apparent increase in variability on the longer timescale may be 
indicative of a physical process like solar activity cycles. 
However, the difference in the two distributions is only margin- 
al and, given the limited number and quality of the measured 
L x values, these results, along with those of Schmitt et al. 
(1993), only hint at the possibility of long-term activity cycles 
for the late-type stars in the Pleiades. 

5. SUMMARY 

Summarizing, we have used three deep ROSAT PSPC 
pointed observations to examine the X-ray spectral and tem- 
poral characteristics of the late-type dwarfs in the Pleiades. We 
have modeled the PSPC pulse-height spectra of individual X- 
ray-bright Pleiades stars with one- and two-temperature RS 
optically thin plasma models. The fits suggest moderately 
hotter coronal temperatures for K and rapidly rotating G stars 
than for F and slowly rotating G stars. For F-, G-, and K-type 
Pleiades stars, we find a general increase in coronal tem- 
perature with Lx/L^,. No M-type Pleiades stars were detected 
with sufficient counts to constrain the spectral models. Our 
findings are in general agreement with those of Schmitt et al. 
(1990) from their analysis of active late-type stars detected by 
the Einstein Observatory . Because the Pleiades contain a 
handful of very rapidly rotating G stars, we have been able to 
show that these G stars have coronal temperatures similar to 
the K stars. 

In order to probe the fainter Pleiades members as well, we 
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have constructed composite X-ray spectra for different spectral 
types as a function of v sin i. Model fits to the composite 
spectra confirm the results based on individual spectral fits. In 
addition, we have determined hardness ratios for nearly all the 
late-type stars. Plots of hardness ratio versus L x and W^bo, 
indicate generally higher ratios (i.e., harder spectra) with 
increased X-ray activity, at least among the F and G stars. 

X-ray light curves have revealed a dozen strong flares in the 
PSPC data. Ten of these flares are associated with known 
Pleiades cluster members. During flares. X-ray luminosities 
rose by factors of 3-20 with rise times in the range 1-10 ks and 
decay times in the range 5-15 ks. The ratio of rise to decay time 
was typically ~ similar to long-duration flares observed by 
EXOSAT and Ginga (Pallavicini et al. 1990; Stern et al. 1992). 
Overall, the Pleiades flares are among the brightest X-ray 
flares ever reported on single, late-type main-sequence dwarfs. 

In their analysis of a single PSPC field of the southwest 
portion of the Pleiades, Micela et al. (1995) detected three 
strong flares. This is consistent with our detection of 12 flares 
in three fields. By considering the number of detected flares, the 
number of observed stars, and the typical monitoring time of a 
single star, we estimate that strong, long-duration (f ~ 10 4 s) 
X-ray flares occur approximately 1 % of the time. This appears 
to be somewhat higher than the flaring frequency observed in 
the Hyades. In fact. Stern et al. (1994) detected only one flare in 
one PSPC field in the Hyades. This may be indicative of the 
overall decrease in activity and rotation rate from the age of 
the Pleiades to the age of the Hyades. On the other hand, 
Pleiades flares occur less often than on the UV Ceti-type flare 
stars (see Pallavicini et al. 1990). 

In fact, Fleming et al. (1995) find that the flaring rate for dMe 
stars in the ROSAT all-sky survey may exceed 10%. Thus, the 
low flaring rate in the Pleiades may appear somewhat sur- 
prising given the youth and high rotational velocities of the 
Pleiades stars. However, the flares we detect in the Pleiades are 
extremely energetic events, with most flare luminosities in 
excess of 3 x 10 30 ergs s -1 . We detect very few less energetic 
flares. Hence, we may be substantially underestimating the 
flaring frequency. 

The duration and peak luminosities of the Pleiades flares are 
reminiscent of strong flares detected on Algol by EXOSAT 
(van den Oord & Mewe 1989) and Ginga (Stern et al. 1992) and 
on T Tauri stars in p Oph (Montmerle et al. 1983) and Orion 
(Gagne et al. 1995). We have modeled one of the flares, associ- 
ated with the late-K Pleiades star H II 1516, with a simple 
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quasi-static cooling loop model. From the PSPC spectrum and 
light curve of the flare, we determine a peak emission measure 
in excess of 10 54 cm -3 and infer a peak electron density n e Z 
10 1 1 cm _ 3 . This suggests a large-scale flare event. 

Since the Pleiades have been extensively monitored by both 
Einstein and ROSAT y we have been able to examine the varia- 
bility of individual stars on timescales of — 1 yr between indi- 
vidual ROSAT exposures and over a decade between the 
Einstein and ROSAT missions. We find that on ^1 yr time- 
scales, approximately one-fourth of the late-type Pleiades stars 
are variable by more than a factor of 2. This is approximately 
the same fraction derived by us for the T Tauri stars in the 
Orion Nebula region (Gagne & Caillault 1994; Gagne et al. 
1995). Based on Einstein observations of the p Oph cloud, 
Montmerle et al. (1983) suggest that the flux variations in T 
Tauri stars may be the result of continual flaring; however, the 
Pleiades ROSAT data are inconclusive. 

On decade-long timescales, up to 40% of the late-type stars 
are variable by more than a factor of 2. This result is in good 
agreement with the fraction determined by Schmitt et al. (1993) 
based on a similar analysis of the ROSA T all-sky survey data. 
Schmitt et al. (1993) have interpreted their findings as possible 
evidence of solar-like activity cycles in the Pleiades. However, 
we warn that the uncertainties involved in the determination of 
count-rate ratios are quite high because of the low detection 
threshold of the X-ray surveys. Moreover, we have compared 
data on only two timescales, and the level of variability is only 
marginally higher on the longer timescales. In order to deter- 
mine the true nature of X-ray variability in late-type stars, 
frequent sampling over a much longer time line is required. 

Note added in manuscript . — A strong flare detected in the 
PSPC data for the Hyades cluster member H II 2411 and 
reported by Preibisch & Neuhauser (1995) was not discussed 
herein because this star is not a Pleiades cluster member. 
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